Some Observations on the TLM Numerical Solution of the Laplace Equation

This paper describes progress on the TLM modelling of the Laplace equation, in particular, how the rate of convergence is influenced by the choice of scattering parameter for a particular discretisation. The hypothesis that optimum convergence is achieved when the real and imaginary parts for the lowest harmonic in a Fourier solution cancel appears to be upheld. The Fourier solution for the problem has been advanced by a better understanding of the nature of the initial excitation. The relationship between the form of the initial condition used in this and many other numerical solutions of the Laplace equation and oscillatory behavior in the results is given a firmer theoretical basis. A correlation between TLM numerical results and those obtained from matrix spectral radius calculations has confirmed much previous work.